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ABSTRACT 

Gyobyu  pipeline,  one  of  the  water  supply  network 
system  of  Yangon  region,  since  this  pipeline  was 
constructed  in  1940,  any  seismic  design  was  not 
considered.  During  the  recent  earthquake,  13.3.2017 
Taikkyi  Earthquake,  Magnitude  5.1,  pipe  supports 
damage  were  observed  at  many  locations  in  the 
Gyopyu  reservoir  area.  Local  site  effect  is  one  of  the 
most  important  aspects  in  the  assessment  of  seismic 
hazard  and  many  of  the  destructive  earthquakes  in  the 
past  have  shown  that  local  site  conditions  have  major 
effects  on  the  ground  shaking.  In  this  study,  the 
recorded  ambient  noise  along  pipeline  was  used  to 
compute  the  H/V  spectral  ratio  and  then  to  evaluate 
the  local  site  effect  of  the  area.  The  HVSR  analysis 
has  been  carried  out  using  Cygwin  Terminal 
Software.  The  predominant  frequency  (F)  ranges 
between  1.17  to  4.41  Hz  and  amplification  factor  (A) 
varies  from  1.07  to  3.25.  The  vulnerability  index  (Kg) 
having  values  ranging  from  0.44  to  7.9.  The 
vulnerability  index  (Kg)  is  indicated  the  soil  damage 
level  due  to  ground  motions. 

KEY  WORDS:  Predominant  frequency,  amplification 
factor,  vulnerability  index,  micro  tremor  HVSR 

I.  INTRODUCTION 

Water  supply  is  vital  to  the  function  of  an 
industrialized  society  and  important  to  emergency 
response  and  recovery  after  earthquakes.  There  are 
four  reservoirs  for  Yangon  water  supply  system.  One 
of  the  reservoirs  is  Gyobyu  reservoir.  The  water 
transmission  pipeline,  however,  in  Yangon,  is 
installed  on  the  ground  supported  by  the  concrete 
foundation  from  the  reservoir  to  the  downtown 
terminal  point  by  segmented  steel  pipe  having 
concrete  lining  with  dresser  coupling  joint.  Gyobyu 
water  pipelines  crossed  six  townships.  They  are 
Taikkyi,  Hmawbi,  Mingalardon,  Mayangon,  Yankin, 
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Figure.  1  Schematic  Configuration  of  Existing  Water 
System  of  Yangon  (JICA,  2012) 

In  the  study,  the  present  water  supply  system  of 
Yangon  which  has  been  operated  over  100  years,  all 
pumping  stations  in  reservoirs  and  main  transmission 
pipes  have  already  aged  and  beyond  its  lifespan, 
because  of  the  long  economic  sanction  to  Myanmar. 
Gyobyu  pipeline  was  constructed  in  1940  and  has 
been  operated  over  70  years.  After  that,  this  pipeline 
suffered  Taikkyi  earthquake  effect  in  March  2017. 
According  to  the  report  of  responsible  person  from 
YCDC,  over  200  nubers  of  pipe  supports  damaged 
between  Gyobyu  reservoir  and  near  Phugyi  reservoir 
(10.5miles).  This  is  the  main  reason  to  select  this  area. 


and  Bahan.  Borehole  data  can’t  get  for  all  townships. 
Therefore,  microtremor  measurement  is  required 
along  the  pipeline. _ 
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Gyobyu  reservoir  is  located  in  Taikkyi  Township.  The 
geological  structure  of  Taikkyi  area  is  underlain  by 
alluvial  deposits  of  Bago  Yoma  fan.  Geographically, 
Taikkyi  is  situated  in  central  Myanmar.  Tectonically, 
it  is  located  in  the  Ayeyarwady  Delta  basin  in  the  west 
of  Sagaing  fault.  Structurally,  the  town  is  located  very 
close  to  a  major  strand  of  the  Sagaing  fault.  [6]  The 
objective  of  this  paper  is  to  investigate  the  soil 
conditions  of  the  selected  area  where  microtremor 
measurements  are  carried  out  along  Gyobyu  lifeline 
within  Taikkyi  area.  The  following  figures  show  the 
pipe  support  cracks  and  their  location  within  Taikkyi 
area.  Each  pipe  independently  moved  vertically 
during  earthquake.  In  recent  earthquake,  any  leakage 
did  not  occur.  The  purpose  of  this  study  is  to  provide 
fundamental  frequency,  amplification  factor  and 
vulnerability  index  values  for  the  study  area. 


Figure.2  Location  of  pipe  supports  damage  within 
Taikkyi 


Figure.  3  Pipe  displacement  after  Taikkyi  earthquake 


Figure.4  some  photos  of  pipe  support  damage 


II.  METHODOLOGY 

In  the  present  study  the  micro  tremor  technique  was 
used  to  illustrate  the  local  site  effect  for  study  area. 
The  single-station  method,  where  a  three-component 
records  form  a  single  seismometer  is  processed  to 
yield  a  spectrum  of  horizontal  to  vertical  spectral  ratio 
and  these  routine  spectral  techniques  can  easily  be 
applied  to  estimate  the  dominant  frequency  of 
vibration  of  the  sedimentary  structure.  These 
frequencies  of  vibration  are  closely  related  to  the 
physical  features  of  the  site  under  study,  i.e.,  layer 
thickness,  densities  and  shear  wave  velocities.  In 
addition,  Microtremor  also  dominated  by  surface 
waves  that  can  be  used  to  determine  the  Rayleigh 
wave  dispersion  curve  without  the  need  for  an 
artificial  source.  Microtremor  analysis  can  be 
performed  using  HVSR  method  (Horizontal  to 
Vertical  Spectral  Ratio).  This  method  is  one  of  the 
easiest  and  most  inexpensive  ways  to  understand  the 
properties  of  the  subsurface  layer  structure  without 
causing  disruption  to  the  structure.  HVSR  method  is  a 
method  that  is  used  as  an  indicator  of  subsurface 
structure,  which  shows  the  relationship  between  the 
ratio  of  the  spectrum  of  H/V  with  the  ellipticity  curve 
of  Rayleigh  waves.  Comparison  of  the  spectral  of  H  / 
V  is  the  ratio  between  the  amplitude  of  Fourier 
spectral  of  horizontal  and  vertical  components  of 
microtremor  [7].  Rayleigh  Ellipticity  is  the  ratio  of  the 
amplitude  of  spectral  components  of  the  horizontal  to 
the  vertical  [8].  Ellipticity  curve  that  is  used  in 
earthquake  engineering  is  to  determine  the 
fundamental  frequency  of  earthquake  locations  and  is 
very  responsive  to  the  location  and  depth  of  a  layer. 
Microtremors  can  be  classified  in  two  major  groups 
according  to  their  frequency  values:  frequency  values 
greater  than  1  Hz  (up  to  20  Hz)  known  as  “short 
period”  which  refer  to  shallow  soil  layers.  The  depths 
can  change  between  0  to  100  meters,  performed  with 
the  effects  of  the  local  surface  sources  such  as 
industry  machines,  wind,  and  traffic  noises. 
Frequencies  lower  than  1  Hz  called  as  “long  period 
microtremor”  which  refers  to  deeper  soil  layers  like 
rigid  rock  formation  and  performed  with  natural 
forces  such  as  ebb  and  tide,  wind  and  underground 
seismic  activities.  Nakamura  (1989)  proposed  a 
method  inferring  site  amplification  factors  to  incident 
seismic  shear  waves  using  microtremor  HV  ratio  at 
single  site.  [9].  In  this  study,  analyzing  the  site 
response  using  the  HVSR  method  and  calculated  Kg 
values  based  on  the  microtremor  data.  Data  was 
reprocessed  using  Cygwin  Terminal  software  to 
estimate  the  Amplification  factors  (A)  and 
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fundamental  frequency  (F)  for  each  observation  point. 
Kg  index  was  estimated  for  each  site  based  on 
Nakamura  (1996,  1997,  2000)  technique.  [10] 


A.  Procedure  for  analysis 

This  is  the  general  procedure  for  calculation  of  H/V 
spectral  ratio  from  record  observation. 


Figure.5  Procedure  of  H/V  spectral  ratio 


After  that,  the  vulnerability  index  values  were 
calculated  in  the  current  study  from  the  microtremor 
measurements.  Seismic  vulnerability  index  (Kg)  is 
one  of  the  simplest  physical  parameter  that  could  be 
used  to  analyze  the  local  geology  condition.  Seismic 
vulnerability  index  was  derived  from  mathematical 
relationship  between  the  resonance  frequency  and 
amplification  factor  of  the  ground  natural  vibration 
that  obtained  in  horizontal  vertical  spectral  ratio 
(HVSR)  curve  analysis.  The  vulnerability  index 
values  are  calculated  from  using  the  equation  (1)  that 
is  expressed  by  Nakamura  (1997): 

Kg=  (1) 

F 

where  Kg  is  the  vulnerability  index,  A  is  the 
amplification  factor  and  F  is  the  fundamental 
frequency. 

Kg  value  reflect  local  site  effect  and  can  be 
considered  as  an  indicator  which  might  be  useful  in 
selecting  weak  point  of  ground  especially  along 
pipeline  areas. 


III.  MICROTREMOR  MEASUREMENTS 
ALONG  PIPELINE 

Microtremor  measurements  were  performed  along 
pipeline  started  in  November  2017  to  still  now.  The 
location  of  all  sites  is  shown  in  Fig.7.  The  Global 
Positioning  System  (GPS)  was  used  to  determine  the 
location  and  elevation  of  the  measurement  sites.  The 
ambient  noise  vibrations  were  recorded  using  LS8800 
equipment.  The  sampling  rate  was  200  Hz  with  a 
recording  duration  of  30  minutes.  The  measurements 
were  taken  during  the  day,  the  distance  between  sites 
is  about  500  meters  and  some  sites  is  about  1 
kilometer  due  to  environmental  constraints.  As  for  the 
measurement  along  pipeline,  over  the  one  hundred 
and  forty  single  station  sites,  which  covered  almost 
the  whole  area  for  pipeline,  were  measured  by  using 
GARMIN  GPS  78  and  LS8800  data  logger  shown  in 
fig-5.  _ 


Figure.6  Observed  microtremor  instrument 


Figure.7  Micro  tremor  measuring  points  along 
Gyobyu  pipeline 
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Figure.  8  Location  of  Microtremor  measuring  points 
within  Taikkyi  area 


beside  pipeline 


A.  H/V  technique  and  H/V  inversion 

The  H/V  analysis  is  based  on  the  theory  first 
introduced  by  Nogoshi  and  Igarashi  [1]  and  further 
developed  by  Nakamura  [2],  who  demonstrated  that 
the  spectral  ratio  between  the  horizontal  and  vertical 
components  of  a  microtremor  record  is  related  to  the 
fundamental  resonant  frequency  at  the  study  site,  the 
application  and  interpretation  of  single  point 
microtremor  measurements  are  very  straight  forward 
and  convenient,  H/V  method  has  become  very  popular 
in  site-specific  investigations  for  study  area. 


In  the  case  of  the  H/V  inversion,  the  trade-off  between 
depths  and  velocities  cannot  be  solved  using  only  the 
H/V  curve  [3].  Therefore,  it  becomes  necessary  to 
retrieve  some  of  these  parameters  (VS  or  thickness) 
from  other  independent  sources.  In  this  work,  reliable 
data  from  borehole  drills  have  been  used  to  define  the 
parameter  space  and  specifically,  to  constrain  the 
number  of  layers  and  their  corresponding  thickness 
ranges.  For  all  the  generated  models,  the  ellipticity  of 
the  fundamental  mode  Rayleigh  wave  is  computed 
(forward  problem)  and  compared  with  the  empirical 
H/V  curve  through  a  misfit  function  [4],  which 
indicates  how  far  the  generated  models  are  from  the 
observed  H/V  results. 


IV.  RESULTS  AND  DISCUSSION 
A.  Distribution  of  soil  natural  frequency  (F)  and 
amplification  factor  (A) 

The  results  of  microtremor  HVRs  for  a  single-station 
measurement  along  pipeline  within  Taikkyi  area 
generally  show  the  following  figure  that  is  one-point 
measurement.  The  fundamental  frequencies, 
amplification  factor  and  vulnerability  index  values  as 
shown  in  Table  I. 


Frequency 


Figure.  10  Example  of  result  of  HVSR  microtremor 
analysis 

And  then,  the  measurement  points  were  shown  in  an 
acceptable  grid  model,  which  was  considered  to 
accomplish  good  quality  contour  maps  (ArcGIS 
software)  that  can  display  the  distribution  of  the 
fundamental  frequency,  amplification  factor  and 
vulnerability.  The  figure.  1 1  shows  the  distribution  of 
frequency  (F)  where  the  distribution  of  natural 
frequency  is  classified  as  equal  interval,  ranging 
between  1.1 7 -4. 41  Hz.  It  could  be  noted  that  frequency 
(F)  is  associated  with  the  depth  of  bedrock.  The 
smaller  of  frequency  (F)  value,  the  greater  of  depth  of 
bedrock. 
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Figure.  1 1  the  distribution  of  fundamental  frequency 
values  in  the  study  area 
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The  figure  12  represents  the  amplification  factor  (A) 
or  peak  ratio  HVSR  spectrum  in  investigation  sites 
ranging  between  1.07  to  3.25.  High  amplification 
factor  (A>3)  is  found  in  some  area.  Different  value  of 
A  might  be  found  in  the  same  value  of  natural 
frequency.  It  can  be  noted  that  the  variation  of  A 
value  is  not  strongly  affected  by  the  soil  depth. 
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Figure.  12  the  distribution  of  amplification  factor 
values  in  the  study  area 
TABLET  FREQUENCY,  AMPLIFICATION 
FACTOR  AND  VULNERABILITY  INDEX  FOR 
STUDY  AREA 


Point 

Nos 

Freque 

ncy 

(F) 

Amplif 

ication 

factor 

(A) 

Vulner 

ability 

index 

(Kg) 

VS 

(30) 

Soil 

Profile 

Type 

1 

2.56 

2.08 

1.69 

354.81 

Sd 

2 

1.61 

1.84 

2.10 

374.26 

Sc 

3 

1.22 

1.85 

2.81 

328.12 

Sd 

4 

1.64 

3.16 

6.09 

267.23 

Sd 

5 

2.78 

3.07 

3.39 

329.44 

SD 

6 

3.71 

2.82 

2.14 

404.25 

Sc 

7 

4.42 

2.69 

1.64 

473.3 

Sc 

8 

2.03 

2.88 

4.09 

355.88 

Sd 

9 

1.51 

2.49 

4.11 

270.6 

Sd 

10 

2.73 

3.19 

3.73 

235.6 

SD 

11 

2.86 

1.93 

1.30 

441.11 

Sc 

12 

2.29 

1.46 

0.93 

515.27 

Sc 

13 

1.88 

1.86 

1.84 

329.47 

Sd 

14 

1.47 

2.25 

3.44 

269.66 

Sd 

15 

2.64 

1.07 

0.43 

333.7 

SD 

16 

2.78 

2.95 

3.13 

397.74 

Sc 

17 

1.68 

2.38 

3.37 

222.99 

SD 

18 

1.24 

1.6 

2.06 

450.5 

Sc 

19 

1.17 

2.72 

6.32 

173.11 

Se 

20 

1.27 

3.17 

7.91 

135.47 

SE 

21 

1.22 

2.96 

7.18 

139.01 

Se 

22 

1.71 

1.87 

2.04 

314.02 

Sd 

23 

1.25 

3.11 

7.74 

234.54 

Sd 

24 

1.19 

2.32 

4.52 

227.27 

SD 

25 

1.56 

2.19 

3.07 

230.77 

Sd 

26 

1.46 

2.05 

2.88 

275.51 

Sd 

27 

1.32 

2.65 

5.32 

169.86 

Se 

28 

2 

2.52 

3.18 

267.61 

Sd 

29 

2.54 

2.31 

2.10 

237.96 

SD 

30 

1.88 

1.73 

1.59 

360 

Sd 

31 

1.66 

3.25 

6.36 

156.13 

SE 

32 

1.46 

2.55 

4.45 

270.61 

Sd 

33 

2.1 

2.69 

3.45 

470.31 

sc 

34 

2.54 

3.02 

3.59 

325.44 

Sd 

B.  Distribution  of  Soil  Vulnerability  Index  (Kg) 

The  use  of  microtremor  was  later  expanded  to  identify 
the  vulnerability  index  for  ground  (Nakamura,  1996, 
1997,  2000,  2009;  Daryono  et  al„  2009)  [11]. 
Nakamura  (1997)  introduced  a  vulnerability  index 
parameter  (Kg),  which  combined  A  and  F  to 
determine  soil  damage  level  due  to  the  local  site 
effects.  Thus,  Kg  can  be  considered  as  an  index  to 
indicated  easiness  of  deformation  of  points  which  is 
expected  useful  to  detect  weak  points  of  the  ground. 
To  estimate  soil  vulnerability  index  (Kg),  the  value  of 
shear  strain  (y)  need  to  be  considered  (Nguyen,  et.al., 
2004).  The  vulnerability  index  (Kg)  is  limited  in  that 
from  y  =  1000  x  10'6  ground  begins  to  show  non¬ 
linear  character  and  in  y  >  10,000  x  10"6  large 
deformation  and  collapse  occur.  Nakamura  (1997) 
had  written  the  equation  2  that  as  follows: 

r=Kgx<Xb  (in  10"6)  &  Kg=A*_(in  10“6)  (2) 

F 

where  ab  is  the  acceleration  of  basement  ground,  A  is 
amplification  factor,  F  is  the  fundamental  frequency 
of  surface  ground.  Based  on  the  above  equation,  the 
effective  strain  can  be  estimated  by  multiplying  Kg 
value  with  maximum  acceleration  of  basement  ground 
ab  in  Gal  (cm/s2)  and  can  be  considered  as  an  index  to 
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indicate  easiness  of  deformation  of  measured  points, 
which  is  expected  to  be  useful  to  detect  weak  points 
of  the  ground.  Figure.  13  shows  the  distribution  of 
vulnerability  index  (Kg)  having  values  ranging  from 
0.44  to  7.9.  In  this  study,  ab  =  0.2  g  (200  Gal)  and  V 
=  800  m/s,  for  y  >1000  x  10"6  then  the  Kg  value  >  5; 
and  for  y  >  10,000  x  10'6  the  Kg  value  >  50.  Based  on 
the  above  criteria,  for  the  value  of  Kg50  was  not 
found  and  Kg  >  5  was  found  at  some  points  in  this 
study  area.  Some  area  are  found  that  vulnerability 
index  (Kg)  <  5.  The  minimum  value  of  (Kg)  is  0.43 
and  maximum  value  of  (Kg)  is  7.9  for  this  study  area. 


Taikkyi 

95°59’0”E  96°0'0”E  96°rO"E  96°2'0"E  96°3’0"E 


Figure.  13  The  distribution  of  vulnerability  index  (Kg) 
values  in  the  study  area 

The  following  figures  show  frequencies,  amplification 
factor  and  vulnerability  (Kg)  value  for  each  site. 
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Fig.  14  (a)  Frequencies  and  Period  (b)  amplification 
factor  (c)  vulnerability  (Kg)  value  for  each  site 


C.  Comparison  between  actual  deformation  and 
measured  results 

In  this  study,  table  1  shows  the  amplification  factor, 
frequency,  shear  waves  and  Soil  Profile  Type  for 
study  area.  It  was  found  that  the  higher  vulnerability 
index  (Kg)  at  the  measured  points  4,  19,  20,  21,  23,  27 
and  31  as  shown  in  Fig.  8.  Almost  of  these  points  are 
soil  profile  type  SE.  The  pipe  supports  damage  were 
found  between  these  points.  The  worst  damages  were 
found  between  5  to  10  miles  far  from  Gyobyu 
reservoir.  The  frequencies  (F)  of  above  6  points  were 
lower  than  1.7  HZ  values  and  the  vulnerability 
indexes  (Kg)  were  great  Kg  >  5.  The  vulnerability 
index  Kg  is  significantly  related  to  the  damage  of  pipe 
supports  when  the  ground  is  not  severely  deformation, 
a  result  of  the  study  of  all  measured  points  are  shown 
in  Fig.  11,  12,  13  and  14. 

V.  CONCLUSION 

The  microtremor  records  along  Gyobyu  pipeline  are 
analyzed  for  34  positions  with  Taikkyi  area.  An 
attempt  has  been  made  to  interpret  variations  of 
microtremor  H/V  spectral  ratios  in  terms  of  the  local 
site  effect.  The  main  conclusions  can  be  summarized 
as  follows: 

>  Seismic  waves  travel  through  different  materials 
at  different  speeds. 

>  Seismic  waves  travel  through  different  materials 
with  different  amplitude. 

>  Bedrock  has  high  frequency  low-amplitude  waves 
(less  destructive). 

>  Sediment  has  low-frequency  high  amplitude  waves 
(more  destructive). 


Therefore,  the  observe  amplitude  and  frequency  are 
reasonable  for  this  concept.  Site  4,  19,  20,  21,  22,  23, 
27,  31  are  more  damage  pipe  support  than  other  site. 
From  the  results,  higher  amplification  is  more 
destructive  at  pipe  support  in  actual  damage 
condition.  There  was  not  found  pipe  support  damage 
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at  that  point.  In  higher  value  of  Kg  (Kg  >  5)  in  some 
points,  there  are  more  damage  of  pipe  supports  other 
than  points  (Kg  <  5)  in  recent  earthquake.  These 
results  show  that  the  weak  location  of  the  soil  for  this 
area,  corresponding  to  recent  earthquake,  can  retrofit 
the  pipe  supports  at  these  location  which  reveal  the 
necessity  for  an  appropriate  emergency  management 
program. 
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